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NOMENCLATURE 
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N 
T 

m,  y(z) 

a. 


ai,  02 
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f 

/c 
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J 

ko,  ki 
m,  p,  q,  r 
s 

J2p+1 

Xo,  Xl,  •  ■  Xff 

yo,  yi . ,  y* 

z 

a 

Oc 


Attenuation  (dB) 

Value  of  A  at  fuft 
(-1  )*C'P 

/^N-r 

Binomial  coefficient,  =Nl/(N—p)\p\ 

£  Mr)BP^r) 

4  —  0 

Numerator  coefficients  of  digital  filter  transfer  function 
Denominator  coefficients  of  digital  filter  transfer  function 
Transfer  function  of  analogue  filter 
Transfer  function  of  digital  filter,  =  Y(z)/X(z ) 

Order  of  filter 
Sampling  period 
2-transform  of  xn,  JMv 

Oc-J»  n  (cos  [(m  -  l)w/2V]/sin  (imr/lN)}  for  p  =  1,2 . N;  (a0  =  1) 

m  —  1 

N  N 

X  Ep  cos  (jxuT),  2  Ep  sin  (ptoT) 

P-0  p- 0 

jy  jy 

2  Fp  cos  (jxoT),  2  Fp  sin  (ptoT) 

p-0  p-0 

Digital  frequency  (Hz) 

Digital  cut-ofif  frequency  (Hz) 

Specified  values  of/  in  passband,  stopband 

<-D,/2 

Constants  used  to  test  for  a  drop-out 
Summation  integers 
Laplace  transform  variable 

Poles  of  H(s)  in  left-half  s  plane  for  p  =  0, 1, ....  N  —  1 

Unfiltered  sequence  of  data 

Filtered  sequence  of  data 

z-transform  variable 

Analogue  frequency  (rad/s) 

Analogue  cut-off  frequency  (rad/s) 


1.  INTRODUCTION 

In  1979  the  Aeronautical  Research  Laboratories  (ARL)  conducted  flight  trials  on  a  Royal 
Australian  Navy  Sea  King  Mk.  SO  helicopter.  The  aim  of  the  trials  was  to  record  data  of  good 
quality  which  could  be  used  to  validate  a  mathematical  model  developed  by  ARL  [1,2].  The 
data  acquisition  system  installed  at  the  time  of  the  trials  allowed  32  channels  of  data  to  be 
recorded  on  magnetic  tape  in  serial  digital  form,  at  sampling  frequencies  of  60,  30,  or  15  Hz  [3]. 
Despite  the  inclusion  of  an  analogue  Butterworth  lowpass  filter  on  most  channels  of  recorded 
data,  it  was  found  that  there  was  significant  noise  remaining  in  a  number  of  these  channels, 
especially  those  measuring  components  of  linear  acceleration  and  angular  velocity.  It  was  decided 
therefore  to  design  and  implement  digital  Butterworth  lowpass  filters  that  could  be  used,  option¬ 
ally,  to  remove  noise  from  any  channel.  This  type  of  filter  was  selected  because  of  its  good  (i.e. 
close  to  ideal)  amplitude  response,  coupled  with  a  reasonably  linear  phase  response,  at  lower 
frequencies  in  the  passband  [4]. 

The  digital  transfer  function  is  usually  represented  either  as  a  series  (i.e.  cascade)  or  parallel, 
or  combination  of  both,  of  second-order  and  first-order  sub-filters  (Ref.  5,  p.  44;  Ref.  6,  p.  26; 
Ref.  7,  p.  43).  Direct  representation  as  a  high-order  filter  is  generally  avoided  for  reasons  of 
numerical  accuracy.  However,  the  direct  form  is  much  simpler  to  program  and  double  precision 
may  be  used  to  retain  sufficient  accuracy.  With  recent  trends  in  computing  towards  efficiency  and 
simplicity  in  programming  effort  at  the  possible  expense  of  CPU  time,  it  was  decided  to  use  the 
direct  form.  Because  no  reference  could  be  found  in  which  the  filter  coefficients  are  derived  for  a 
general-order  filter  in  the  direct  form,  the  lengthy  algebraic  process  is  included  in  Section  2  in  a 
description  of  the  basic  theory  used  to  design  the  digital  filter. 

In  Section  3,  flight  data  recorded  for  a  particular  manoeuvre,  together  with  power  spectra, 
are  presented  for  most  channels  that  were  properly  functioning  at  the  time  of  recording.  In 
obtaining  power  spectra,  the  low  frequency  signal  components  are  first  removed  to  accentuate 
unwanted  frequency  components  at  and  above  the  one-per-rev  rotor  frequency  of  3  •  5  Hz. 
In  Section  4,  the  filtering  procedure  adopted  for  each  channel  is  given.  In  cases  where  the  control 
input  and  responses  to  it  have  a  substantial  energy  content  at  relatively  high  frequencies,  the 
analogue  filtering  may  cause  some  distortion  of  the  wanted  signal.  Hence,  care  should  be  taken 
to  avoid  further  distortion  when  using  a  digital  filter.  For  this  reason  and  because  of  the  absence 
of  significant  noise,  control  inputs  (i.e.  pilot  stick  and  pedal  displacements)  and  displacement 
quantities  directly  linked  (i.e.  mechanically  or  electrically  rather  than  aerodynamically)  to  these 
inputs  are  not  filtered.  For  all  channels  requiring  filtering,  two  standard  fifth-order  filter  speci¬ 
fications  are  given.  One  has  a  high  attenuation  of  50  dB  at  the  rotor  frequency  of  3  •  5  Hz,  giving 
a  low  cut-off  frequency  of  1  •  12  Hz,  and  removes  noise  at  and  above  the  rotor  frequency.  It  is 
suitable  for  variables  where  the  energy  is  confined  to  relatively  low  frequencies,  and  its  effective¬ 
ness  is  demonstrated  on  three  different  types  of  signal.  The  other  filter,  which  has  a  higher  cut-off 
frequency  of  4  Hz,  is  used  alternatively  where  the  high  frequency  content  is  significant.  Generally, 
of  the  kinematic  quantities  measured,  the  only  one  where  this  is  likely  to  be  the  case  is  vertical 
acceleration  when  there  is  a  violent  manoeuvre.  Hence,  for  such  a  manoeuvre,  the  effect  of  both 
filters  on  the  vertical  acceleration  signal  is  shown. 

2.  DIGITAL  FILTER  DESIGN 

2.1  General  Realization  of  a  Recursive  Filter 

For  a  recursive  realization  of  a  digital  filter  transfer  function,  the  functional  relationship 
between  the  input  sequence  xo,  xi, . . .,  xn  and  the  resulting  output  sequence  yo,  yi,  . .  ys  is 
given  by  the  difference  equation  (Ref.  6,  p.  22;  Ref.  7,  p.  40;  Ref.  8,  p.  197) 

N  N 

£  fpTjr-p  ®  ^.FpyN-p  (1) 

1-0  p=0 
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1 


where  Eo,  Eu  . . E\  and  Fo,  Fi . are  the  filter  coefficients.  Rewriting  Equation  (1)  to 

solve  for  yx  explicitly, 


1  f  v  Jv  “I 

yx  =  7T  Z  EpXN-p  —  Z  Fpyx-p  L  (2) 

/■oLp  =  o  p=  i  J 

The  digital  filter  is  then  realized  in  terms  of  the  constant  filter  coefficients,  which  must  be 
determined. 

In  digital  filter  theory,  it  is  convenient  to  transform  the  above  time-domain  equation  to  the 
z-domain  using  the  z-transform,  in  which  each  uniform  time  delay  in  both  input  and  output 
sequences  is  represented  by  a  multiplication  by  z~l  (Ref.  8,  p.  45).  Hence 

Y(z)  =  Z  Epz~PX(z)  -  I  Fpz-P  T(z)l  (3) 

/'oLp  =  o  p=t  J 

where  X(z)  and  Y(z)  are  the  z-transforms  of  xn  and  yn  respectively.  From  Equation  (3),  the 
transfer  function  from  which  the  digital  filter  may  be  realized  is  given  by 

N 


H(z)  = 


m 

X(z) 


Z  EpZ-r 

P  =  0 
N 


HFpZ-P 
P  =  0 


(4) 


2.2  Butterworth  Digital  Filter  Coefficients 

The  squared  magnitude  response  for  an  analogue  Butterworth  lowpass  filter  is  of  the  form 
(Ref.  8,  p.  211) 

l"0fl)i2=;-7T7 A.^.v  (5) 


1  +  (mile)™ 

whfre  £2  is  analogue  frequency  in  ra d/s,  flc  is  the  cut-off  frequency  (i.e.  amplitude  response  is 
3  dB  down  at  £1  —  ilc),  and  N  is  the  order  of  the  filter.  Section  2.3  shows  how  both  Qc  and  N, 
which  are  assumed  known  in  this  section,  may  be  determined  from  a  design  specification. 

On  substituting  s  =  j£l  in  Equation  (5), 

(6) 


1  +  (j//nc)2" 

The  poles  of  Equation  (6)  are  then  at  equally  spaced  points  on  a  circle  of  radius  Oc  in  the  j-plane. 
By  associating  the  poles  in  the  left-half  j-plane  with  H(s), 


m  ri(s  -  s2p+1) 

p-0 


(7) 


where  (Ref.  9,  p.  494;  Ref.  10,  p.  183) 

S2J)+1  =  Cl^ZP^XMiN  p  =  0,  1,  . 

Alternatively,  Equation  (7)  may  be  written  as  (Ref.  9,  p.  494) 

1 


N  —  1 


where 


H(s)  =  " 

T.aj>sp 

p  =»0 


(8) 


_  P  fcos  [(m  —  1  )rr/2iV 

Op  =  Qc-P  n  \ - “  p  =  1,  2 . N 

m  =  i(  sin  {mn/2N)  j 


ao  =  1 

To  form  the  digital  transfer  function  H(z),  the  bilinear  transformation  is  applied  in  which 
(Ref.  8,  p.  207) 

2(1  -z->) 


s  = 


7X1  +  z-i) 


(9) 


is  substituted  in  Equation  (8)  to  give 


—  Z~ 


where  T  if  the  sampling  period.  Expansion  of  Equation  (10)  in  the  direct  form  of  Equation  (4) 
is  generally  avoided  because  the  effect  of  round-off  errors  in  the  filter  coefficients  becomes 
increasingly  significant  for  higher-order  filters  (Ref.  7,  p.  41).  Therefore,  the  usual  method  for 
expanding  Equation  (10)  is  to  write  the  expression  in  either  series  (cascade)  or  parallel  form  in 
terms  of  second-order  and  first-order  sections  (Ref.  5,  p.  44;  Ref.  6,  p.  26;  Ref.  7,  p.  43).  However, 
the  direct  form  of  Equation  (4)  is  much  simpler  to  program  and  the  problem  of  coefficient 
sensitivity  can  be  overcome  by  using  double  precision  arithmetic  in  the  computer  program  where 
appropriate.  For  overall  efficiency,  the  direct  form  was  chosen,  and  because  no  reference  could 
readily  be  found  in  which  Equation  (10)  is  expanded  in  direct  form  for  a  general  value  of  N,  the 
lengthy  algebraic  process  is  now  given. 

Rewriting  Equation  (10), 

(1  +  rW 

H(z)  =  -N - -  (11) 

ZMz/TW  -  z- i)*(l  +  2-‘)*-p} 
p- o 

Using  the  Binomial  theorem  to  expand  powers  of  (1  +  z_1)  and  (1  —  z~x). 


£  M2/7T I  {( ~l)»z-*C'p]  E 


where 


CN  -  •  m 
’  ( N-pV.pl 


0  N 


=  0  otherwise 

It  should  be  noted  that  the  summation  variables  in  the  denominator  of  Equation  (12)  are  re¬ 
arranged  at  this  stage  in  order  to  arrive  at  the  final  form  of  H{z)  with  the  desired  summation 
variable  p. 

Defining 

MO  =  (-l)pC£ 

Mr)  -  C 

then  in  the  denominator  of  Equation  (12), 

£  l(-l)Pz-PC'J'j:iz-<C?-']  =  ilz-’AMVZlz-'BAr)] 

J>  •»  O  q  =  0  p  -0 

=  [7-°y<o(r)  +  z-'Ai(r)  +  . . .  +  z^AM^B^r)  +  z-'Bi(r)  +  . . .  +  z^^Bk-At)] 

=  z-°[/fo(r)Ro(r)]  +  z~x[A^r)Bi(r)  +  Ai(r)B*r))  +  ...+  z^[A^r)BN(r)  +  Ai[r)Bs-i(r) 
+  • . .  +  Ax(r)3o(r)] 

=  I  z-’Djtr) 

p- o 

where 

Dp(r)  =  E  Aq(r)Bp-q(r) 

1  =  0 

In  the  above  expansion,  terms  are  included  for  convenience  that  may  be  zero,  e.g.  BN(r)  =  0 
if  r  >  0. 

The  denominator  of  Equation  (12)  is  then  given  by' 

£{aA2ITy£z~PDp(r))  =  a^l/Tnz^D^O)  +  z-‘D,(0)  +  . . .  +  z-*D*(0)] 

f  —  0  r-0 


+  ai(2T)i[z-0A*l)  f  z-'Z)i(l)  +  . . .  +  z^DMO] 

+  . .  — <3/7"  °A>(/V)  +  z~xDi(N)  +  . . .  +  z-*Dn{N)] 


=  ^{z~^MVTYDp(r)} 

p  =  0  r  =  0 

Hence,  Equation  (12)  may  be  written  as 


H(z)  =  iv- 


N 

£z-rC* 
p  =  o 


/» 


£  [rPS^/rrJ 

p=0t  r=0  1 

which  is  the  direct  form  of  Equation  (4)  with  the  filter  coefficients  given  by 


(13) 


Ep  =  CNp 

Fp  =  I  a,(2jT)T  Dp{r) 

r  =  o 

Fortran  subroutines  for  obtaining  these  coefficients  are  listed  in  Appendix  A.  Only  the  de¬ 
nominator  coefficients  need  be  calculated  using  double  precision,  since  the  numerator  coefficients 
are  integer. 


2.3  Design  Specification 

For  a  digital  Butterworth  lowpass  filter,  defined  in  the  form  of  Equation  (13),  values  of  the 
parameters  Oc  and  N  that  enable  the  desired  specifications  to  be  met  need  to  be  determined. 
Frequently,  these  design  specifications  are  given  in  terms  of  a  maximum  attenuation  A\  (dB) 
at  a  frequency  f\  (Hz)  in  the  passband  and  a  minimum  attenuation  Az  (dB)  at  a  frequency  fz  (Hz) 
in  the  stopband  (Ref.  8,  p.  214).  To  correct  for  distortion  inherent  in  using  the  bilinear  transfor¬ 
mation,  the  specified  digital  frequencies  f\  and  /a  in  Hz,  after  first  converting  to  radian  frequencies 
<oi  and  <02,  must  be  prewarped  to  the  corresponding  analogue  frequencies  tii  and  ft2  in  rad/s. 
For  general  values  /,  <o,  and  Q,  this  is  achieved  using  the  relations  (Ref.  8,  p.  208) 


<o  =  2v f  (14) 


From  the  definition  of  attenuation  A  in  dB, 

A  =  —20  logio  \H(jCi)\  (16) 

Substituting  for  H(JQ)  in  Equation  (16)  from  Equation  (5), 

Qc  =  £2(10x/ 10  -  l)"1/ZAr  (17) 

On  using  Equation  (17)  at  the  two  specified  frequency  values,  £2C  may  be  eliminated  from  the 
resulting  two  simultaneous  equations  to  give 

logio  [(10^‘°  -  1)/(10^*°  -  1)] 

2  logio  (Qi/Oa)  1  J 

The  value  of  N  given  by  Equation  (18)  is  then  rounded  up  to  the  nearest  integer  and  Oc  is  obtained 
from  either  simultaneous  equation  using  this  value  of  N. 

Unless  the  value  of  N  obtained  from  Equation  (18)  is  already  an  integer,  only  one  of  the 
design  points  may  be  matched  precisely.  The  particular  simultaneous  equation  used  depends  on 
which  of  the  design  specifications  one  wishes  to  match.  It  has  been  found  useful  to  be  able  to 
choose  either  of  these  points  as  a  means  of  identifying  the  filter  with  its  main  purpose.  For 
example,  one  may  wish  to  identify  the  passband  frequency  as  the  cut-off  frequency  (i.e./i  =/c 
at  Ai  =  3  dB),  or  the  stopband  frequency  as  the  rotor  frequency  (e.g.  A2  =  50  dB  at/2  =  3-5 
Hz).  Therefore,  the  above  procedure  is  used  only  to  obtain  an  initial  estimate  for  N.  The  filter 
is  then  specified  in  terms  of  a  selected  value  of  N  and  an  attenuation  at  a  specified  frequency 
anywhere  in  the  frequency  domain.  Equation  (17)  then  yields  £!c  directly. 


r 


r  H 
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2.4  Phase  Delay 


Ideally,  it  would  be  desirable  to  completely  eliminate  the  non-linear  phase  response  of  the 
Butterworth  filter  from  the  filtered  data.  One  way  of  achieving  this  is  to  process  the  data  in  both  a 
forward  and  backward  direction  using  the  same  digital  filter  (Ref.  11,  p.  194),  which  results  in 
exact  cancelling  of  the  two  phase  shifts.  When  using  this  technique,  allowance  should  be  made 
for  the  squaring  of  the  transfer  function.  Preliminary  investigations  of  data  to  be  filtered  showed 
that  this  added  complication  was  not  warranted,  since  adequate  correction  could  be  made  by 
assuming  linear  phase  shift  and  hence  constant  phase,  or  time,  delay.  Given  the  transfer  func¬ 
tion  H(z),  an  expression  giving  this  correction  is  now  derived. 

In  the  analogue  filter,  frequency  is  measured  along  the  imaginary  axis  in  the  complex 
5-plane,  whereas,  in  the  digital  filter,  frequency  is  measured  along  the  circumference  of  the  unit 
circle  in  the  z-plane  (Ref.  5,  p.  43).  Hence,  on  substituing  z  =  e^T  in  Equation  (4), 

N 


H(e^T)  = 


£  Epe-P*“T 
p  =  o _ 

N 

£  Fve-P^T 
f=o 


ai  —  jaz 
bi  —  jbi 


(19) 


where 

N  N 

ai  =  X  Ev  cos  (P^T)  a2  =  £  Ep  sin  (pwT) 

/>“  o  P  =  o 

N  N 

bi  =  X  Fp  c°s  (/*« 7")  bi  =  £  Fp  sin  ( pwT ) 

/>-  o  />=  o 

Both  the  amplitude  and  phase  response  are  readily  obtained  from  Equation  (19),  with  the  phase 
response  <f>(e>uT)  given  by 

(oibz  —  o«bi\ 
aibi  +  azbz) 

The  phase  delay  t(^<"t)  is  then,  in  units  of  seconds  (Ref.  12,  p.  66) 

r(e'"T)  =  -<fte)“T)/w  (21) 


<t>(eiuT)  =  arctan  | 


(20) 


For  the  digital  Butterworth  lowpass  filter,  the  phase  response  is  reasonably  linear  at  lower 
frequencies  in  the  passband.  Hence,  the  phase  delay,  which  is  assumed  constant,  is  evaluated 
at  =  <i>c/10. 


3.  SPECTRAL  ANALYSIS  OF  FLIGHT  DATA 

The  data  acquisition  system  (DAS)  used  for  the  flight  trials  [3]  allowed,  without  sub¬ 
multiplexing,  20  channels  of  data,  with  12-bit  word  size,  to  be  recorded  at  a  sampling  frequency 
of  60  Hz.  However,  by  sub-multiplexing  each  of  three  channels  to  record  two  separate  quan¬ 
tities  at  30  Hz,  and  each  of  another  three  to  record  four  separate  quantities  at  1 5  Hz,  a  total  of 
32  channels  was  actually  recorded.  Signals  are  recorded  from  the  channels  in  a  regular  sampling 
order  at  equi-spaced  time  increments  throughout  each  complete  sampling  time  interval  of 
l/60th  second.  However,  it  is  assumed  that  all  signals  recorded  within  this  sampling  interval  are 
at  the  clock  time  recorded  at  the  beginning  of  the  interval.  Hence,  a  maximum  relative  time 
delay  error  of  close  to  l/60th  second  is  tolerated.  For  the  sub-multiplexed  channels,  the  particular 
channel  being  recorded  at  each  complete  sampling  interval  is  determined  by  the  two  least  sig¬ 
nificant  bits  of  the  clock  time  word  [3].  Values  for  channels  not  recorded  at  a  specific  time, 
because  of  sub-multiplexing,  are  obtained  by  linear  interpolation,  which  allows  a  common 
sampling  frequency  of  60  Hz  to  be  assumed  when  filtering  the  data. 

Having  converted  the  data  to  a  format  suitable  for  input  to  a  DEC  PDP-IO  computer,  the 
first  stage  in  the  data  refinement  process  is  to  correct  for  drop-outs,  generally  caused  by  a  mis¬ 
placed  significant  bit  in  a  12-bit  binary  word.  For  an  unfiltered  sequence  of  signals  x0,  *i,  xz, 
which  includes  interpolated  values  for  sub-multiplexed  channels,  a  drop-out  is  considered  to 
have  occurred  when 

\X2  —  *l|  >  ko  +  fcll-Xl  -  *ol 

where  k0  and  ki  are  constants  determined  appropriately  by  inspection  of  the  flight  data.  For  a 
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12-bit  word,  giving  an  integer  range  of  0  to  409S,  values  found  appropriate  were  k0  =  100  and 
ki  =  5.  When  a  drop-out  occurs,  the  current  value  is  reset  to  its  previous  values,  i.e.  xz  =  xi. 

The  next  stage  is  to  correct  for  phase  shift  introduced  by  the  tie  of  analogue  Butterworth 
lowpass  filters  in  the  DAS.  Sixth-order  filters  were  used,  with  cut-off  frequencies  of  3,  6,  and  12 
Hz,  generally  corresponding  to  the  sampling  frequencies  of  15,  30,  and  60  Hz.  A  description  of 
the  quantities  measured,  together  with  their  sampling  frequency  and  analogue  filter  cut-off 
frequency,  are  given  in  Table  1.  For  the  purpose  of  obtaining  data  suitable  for  validation  of  the 


TABLE  1 

Quantities  Measured  by  Data  Acquisition  System 


Channel 

number 

Quantity 

measured 

. 

Sampling 

frequency 

(Hz) 

Analogue 

filter 

cut-off 

frequency 

(Hz) 

1 

Cyclic  stick  position — pitch 

60 

12 

2 

Cyclic  stick  position — roll 

60 

12 

3 

Collective  stick  position 

60 

12 

4 

Angle  of  att  ack 

60 

12 

5 

Fore-aft  push-pull  rod  position 

60 

3t 

6 

Lateral  push-pull  rod  position 

60 

12 

7 

Collective  push-pull  rod  position 

60 

12 

8 

Pitch  rate 

60 

12 

9 

Sideslip  angle 

60 

12 

10 

Roll  attitude 

60 

12 

11 

Roll  rate 

60 

12 

12 

Longitudinal  acceleration 

60 

12 

13 

Lateral  acceleration 

60 

12 

14 

Vertical  acceleration 

60 

12 

15 

Pitch  attitude 

30 

6 

16 

Yaw  pedal  position 

30 

6 

17 

Yaw  rate 

30 

6 

18* 

Yaw  attitude 

30 

— 

19 

Lateral  cable  angle 

30 

6 

20 

Longitudinal  cable  angle 

30 

6 

21 

Longitudinal  velocity 

15 

6 

22 

Lateral  velocity 

15 

6 

23 

Dynamic  pressure 

15 

12 

24 

Radio  altitude  (raw) 

15 

3 

25 

Radio  altitude  (smooth) 

15 

3 

26 

Absolute  pressure 

15 

3 

27 

Yaw  push-pull  rod  position 

15 

12t 

28 

Ambient  temperature 

15 

— 

29 

Torque — Engine  1 

15 

— 

30 

Rotor  r.p.m. 

15 

— 

31 

Towed  probe  dynamic  pressure 

15 

3 

32 

Towed  probe  differential  pressure 

15 

3 

*  Switch  selectable  alternative  measurement  of  “Torque— Engine  2". 
t  Filters  inadvertently  interchanged. 

Note:  Acceleration  measurements  include  gravitational  effects. 


Sea  King  mathematical  model  [1,2],  accurate  correction  of  phase  shift  is  considered  necessary 
only  at  lower  frequencies  in  the  passband  of  each  filter.  In  this  region,  the  filters  used  have  a 
reasonably  linear  phase  response  (Ref.  12,  p.  1 12)  so  that  a  constant  time  delay  may  be  assumed. 
Consistent  with  inaccuracies  already  accepted  for  relative  time  delays,  these  corrections  are 
expressed  in  terms  of  complete  sampling  time  intervals.  For  the  filters  with  cut-off  frequencies 
3,  6,  and  12  Hz,  the  corresponding  time  delays  are  12,  6,  and  3  sampling  intervals. 

Flight  data  obtained  for  a  particular  manoeuvre  (ARL  Sea  King  Flight  Record  No.  19100) 
are  presented  in  Figures  1  to  3,  in  which  the  data  is  corrected  for  analogue  filter  phase  delay 
and  converted  to  suitable  engineering  units  using  appropriate  calibration  factors  and  offsets. 
The  channels  are  grouped  into  three  types,  which  include: 

(1)  control  inputs  and  directly  linked  displacements; 

(2)  angular  displacements  and  velocities;  and 

(3)  linear  velocities  and  accelerations,  and  miscellaneous  quantities. 

Because  of  a  temporary  fault  in  the  recording  system,  sideslip  angle  was  not  measured  for  this 
particular  flight,  but  in  other  flights,  the  characteristics  of  the  signal  were  found  to  be  similar 
to  those  of  the  angle  of  attack  signal.  The  incremental  plotter,  used  to  present  the  data  in  Figures 
1  to  3,  plots  in  increments  of  0  01  inch  in  both  X  and  Y  directions.  Hence  to  avoid  distortion  of 
the  signal  along  the  time  scale  caused  by  rounding  equi-spaced  intervals  of  l/60th  second  to  the 
nearest  plot  increment,  the  scale  is  chosen  such  that  each  sampling  interval  is  an  integer  multiple 
of  the  plot  increment  (in  this  case,  the  integer  multiple  is  one,  i.e.  100  samples  per  inch).  Though 
the  duration  of  the  recorded  flight  data  was  18-667  s,  only  the  first  8-333  s  (i.e.  500  sampling 
intervals)  are  shown. 

In  order  to  determine  the  relative  energy  levels  in  the  signals  for  frequencies  up  to  the 
Nyquist  frequency,  which  is  half  the  sampling  frequency,  spectral  analysis  is  used  [13].  In  using 
the  Fast  Fourier  Transform  (FFT)  method  to  obtain  the  Fourier  coefficients  [14],  the  number 
of  sample  values  is  required,  by  convention,  to  be  a  power  of  two.  The  quantity  of  data  available 
for  the  flight  (1121  samples),  allows  2 10  (=1024)  samples  to  be  analyzed  when  the  sampling 
frequency  is  60  Hz.  Since,  for  sub-multiplexed  channels,  interpolated  values  are  not  true  sample 
measurements,  they  should  not  be  used.  Hence,  the  corresponding  number  of  samples  at  30 
and  15  Hz  are  2®  (=512)  and  28  (=256).  The  signal  components  at  low  frequencies  up  to  about 
1  Hz  are  first  removed  by  subtracting,  from  the  unfiltered  data,  values  obtained  using  the  low 
cut-off  frequency  filter  eventually  adopted  (see  Section  4).  This  accentuates  unwanted  frequency 
components  at  and  above  the  rotor  frequency  of  3  •  5  Hz.  The  noise  power  spectrum  is  then 
obtained  by  first  smoothing  with  a  Hanning  window  and  then  smoothing  further  with  weights 
0- 16,  0-68,  and  0- 16  as  defined  by  Blackman  (Ref.  15,  p.  136). 

For  the  control  inputs  and  directly  linked  displacements  in  Figure  4,  there  is  an  absence  of 
significant  noise.  In  Figures  5  and  6,  significant  noise  is  indicated  for 

(a)  yaw  attitude  at  3-5  Hz; 

( b )  roll  rate  and  all  three  accelerations*  at  3  ■  5  Hz  and  higher  harmonics  (including  fifth, 
which  is  blade  passing  frequency  of  17-5  Hz); 

(c)  pitch  rate  at  17-5  Hz; 

( d )  angle  of  attack,  dynamic  pressure,  and  torque  at  unspecified  frequencies  (i.e.  random 
noise);  and 

(e)  rotor  r.p.m.  at  3  -5  Hz  and  higher  frequencies. 

Pitch  attitude,  roll  attitude,  yaw  rate,  longitudinal  velocity,  and  lateral  velocity  do  not  have 
significant  noise. 

*  Though  the  energy  level  for  vertical  acceleration  is  very  high  at  17-5  Hz,  there  is  still 
significant  noise  at  3  •  5  Hz. 


Repetition  of  the  spectral  analysis  procedure  on  other  flight  manoeuvres  produced  very 
similar  results.  Hence  it  was  considered  appropriate  to  draw  general  conclusions  on  filtering  of 
each  channel  from  the  trial  analysed  here. 

4.  FILTERING  OF  FLIGHT  DATA 

In  adopting  a  filtering  procedure  for  each  channel,  one  must  first  decide  on  the  relative 
merits  of  filtering  the  channel.  In  cases  where  the  high  frequency  energy  content  is  significant, 
filtering  is  likely  to  distort  the  signal  to  some  degree.  This  is  because  the  filtering  does  not  distin¬ 
guish  between  wanted  and  unwanted  signal  components. 

Control  inputs  and  directly  linked  displacements  frequently  contain  substantial  high 
frequency  energy,  but  the  noise  component  is  generally  insignificant.  Hence,  these  channels  are 
not  filtered.  Except  in  special  cases,  such  as  violent  manoeuvres,  the  high  frequency  energy 
content  of  other  channels  primarily  arises  from  noise.  Since  most  noise  is  at  and  above  the  rotor 
frequency  of  3  •  5  Hz,  it  was  decided  therefore  to  filter  all  these  channels,  including  those  with 
little  noise,  using  a  fifth-order  filter  with  an  attenuation  of  50  dB  at  3-5  Hz,  giving  a  low  cut¬ 
off  frequency  of  112  Hz.  The  effect  of  the  filter  on  three  types  of  signal  is  shown  in  Figure  7, 
where,  for  clarity,  the  filtered  signals  are  offset  from  the  unfiltered  ones.  It  can  be  seen  that  the 
filter  copes  well  with  the  first  two  types  shown,  where  noise  is  clearly  identifiable  at  specific 
frequencies,  as  well  as  the  third  type,  which  has  random  noise. 

Examples  of  violent  manoeuvres  include  engine  cuts  and  sudden  large  movements  in  the 
collective  stick,  which  generally  result  in  sudden  large  changes  in  vertical  acceleration,  altitude, 
torque,  and  rotor  r.p.m.  For  such  cases,  when  filtering  is  necessary,  a  fifth-order  filter  is  specified 
with  a  higher  cut-off  frequency  of  4  Hz.  The  effect  of  both  filters  on  the  vertical  acceleration  signal 
is  shown  in  Figure  8  for  a  typical  violent  manoeuvre,  with  filtered  values  again  offset  from  un¬ 
filtered  ones.  The  low  cut-off  frequency  filter  clearly  distorts  the  acceleration  response,  while  the 
high  cut-off  frequency  filter  does  not  distort  the  signal  perceptibly.  However,  for  the  latter,  while 
the  main  noise  component  at  17  5  Hz  is  effectively  removed,  significant  noise  is  still  present  at 
3  -5  Hz.  The  transfer  function  coefficients  and  phase  delay  of  each  filter  are  given  in  Table  2. 


TABLE  2 

Specifications  of  Standard  Digital  Butterworth  Lowpass  Filters 


Quantity 

Filter  1 

(low  cut-off  frequency) 

Filter  2 

(high  cut-off  frequency) 

T  (s) 

1/60 

1/60 

fi 

5 

5 

/(Hz) 

3-5 

4-0 

A  (dB) 

500 

3-0 

fc  (Hz) 

112 

4-0 

Qc  (rad/s) 

7  •  033 1 292892882D +00 

2-55189033061 18D+01 

r(s) 

0-459 

0-127 

£o ,  F0 

1,  1 -7478224133541 D+06 

1,  4-5580802280576D+03 

EuFi 

5,  —  8  ■  07698 1 4534886D + 06 

5,  - 1 -662 1 49955 1730D+ 04 

e2,  f2 

10,  1  •  4953389527402 D+ 07 

10,  2  ■  487762821 8574  D+ 04 

£3 » Ea 

10,  —  1 -3862462739192D+07 

10,  —  1 -899I475064559D+ 04 

£4,  £4 

5,  6  •  434534069538 3  D + 06 

5,  7-3691 389324960D + 03 

£5,  £s 

1,  —  1  •  1 962698 1 76 138D+ 06 

, 

1,  —  1  ■  1598727628397D+03 

5.  CONCLUDING  REMARKS 

For  the  purpose  of  removing  significant  noise  from  flight  data  on  a  Sea  King  Mk.  50  heli¬ 
copter,  suitably  designed  digital  Butterworth  lowpass  filters  are  demonstrated  to  be  effective. 
For  simplicity  and  economy  of  programming  effort,  the  direct  form  of  representation  of  the 
transfer  function  as  a  high-order  filter  is  successfully  implemented,  using  double  precision  to 
avoid  the  usual  problem  of  numerical  accuracy  caused  by  high  coefficient  sensitivity.  That  no 
reference  in  the  literature  could  readily  be  found  in  which  the  filter  coefficients  are  derived  in 
direct  form  for  a  general-order  filter  is  probably  because  of  the  generally  accepted  practice  of 
representing  a  high-order  filter  appropriately  as  a  number  of  second-order  and  first-order  sub¬ 
filters.  Hence,  the  lengthy  algebraic  process  is  included  here.  Using  representative  flight  data, 
spectral  analysis  is  used  to  clearly  identify  unwanted  frequency  components.  To  accentuate  these 
components,  which  are  mostly  at  the  rotor  frequency  of  3  -5  Hz  and  higher  harmonics,  the 
signal  components  at  low  frequencies  up  to  about  1  Hz  are  first  removed  by  subtracting,  from  the 
unfiltered  data,  values  obtained  using  the  low  cut-off  frequency  filter  eventually  adopted. 
Channels  measuring  control  inputs  and  directly  linked  displacement  quantities  are  not  filtered 
because  of  the  absence  of  significant  noise  and  to  avoid  further  distorting  the  signal.  Two 
standard  fifth-order  filters  are  shown  to  be  adequate  for  all  filtering  requirements.  One,  with  a 
low  cut-off  frequency  of  M2  Hz,  is  suitable  for  variables  where  the  energy  is  confined  to 
relatively  low  frequencies,  while  the  other,  with  a  higher  cut-off  frequency  of  4  Hz,  is  used  alterna¬ 
tively  where  the  high  frequency  content  is  significant.  For  both  filters,  adequate  correction  for 
phase  shift  is  obtained  on  assuming  a  constant  time  delay. 
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APPENDIX  A 


Fortran  Subprogram  Giving  Transfer  Function  Coefficients  in 
Direct  Form  for  a  Digital  Lowpass  Butterworth  Filter 


SUBROUTINE  TRANSF ( N , SFREQ , OMEGAC ,E,F, A,B) 

C 

C  CALCULATES  NUMERATOR  AND  DENOMINATOR  COEFFS  OF  TRANSFER  FUNCTION 
C  (I.E.  ARRAYS  'E'  AND  'F')  THAT  DEFINE  A  DIGITAL  LOWPASS  BUTTERWORTH  FILTER 
C  WITH  »N'  POLES,  ANALOGUE  CUT-OFF  FREQ  'OMEGAC',  AND  SAMPLING  FREQ  'SFREQ'. 
C  ARRAYS  'A'  AND  'B',  THOUGH  USED  ONLY  LOCALLY  WITHIN  SUBN,  ARE  DIMENSIONED 
C  EXTERNALLY.  THE  TRANSFER  FUNCTION  IS  IN  THE  DIRECT  FORM 
C 

PARAMETER  PID2=1 .570796326794897D0  !  PI/2 

INTEGER  SFREQ  !  1/T 

INTEGER  BCOEFF,E,A,B,D 

INTEGER  P,Q,R,PP,QQ,RR 

DOUBLE  PRECISION  OMEGAC, F, PROD, AR 

DIMENSION  E(0:N) ,F(0:N) ,A(0:N) ,B(0:N) 

C 

DO  210  P=0,N 
PP=P 

E(P)=BCOEFF(N,PP)  !  NUMERATOR  COEFF 

C 

F(P)=0D0 

C 

DO  200  R=0,N 
RR=R 

PROD=1DO 

IF(R.EQ.O)  GO  TO  120 
C 

DO  110  M=1 ,R 

PROD=PROD*DCOS ( (M-1 )»PID2/N)/DSIN(M»PID2/N) 

1 10  CONTINUE 

120  AR=PROD/OMEGAC**R 

C 

DO  130  Q=0,N 

A(Q)=0 

B(Q)=0 

130  CONTINUE 

C 

DO  140  Q=0,R 
QQ=Q 

A(Q)=BC0EFF(RR,QQ)*(-1 )”Q 
140  CONTINUE 

C 

DO  150  Q=0 , N-R 

QQ=Q 

B(Q)=BCOEFF(N-R,QQ) 

150  CONTINUE 

C 

D=0 

DO  160  Q=0,P 
D=D+A(Q)*B(P-Q) 

160  CONTINUE 

C 

F(P)sF(P)+AR»D»(2«SFREQ)*»R  !  DENOMINATOR  COEFF 
200  CONTINUE 

C 

210  CONTINUE 

C 

RETURN 

END 


o  o  n  o  o  non 


APPENDIX  A  (cont) 


INTEGER  FUNCTION  BCOEFF(N,K) 

GIVES  BINOMIAL  COEFFS 
INTEGER  FACT 

BCOEFF =FACT (N)/( FACT (K)*FACT(N-K) ) 

RETURN 

END 


INTEGER  FUNCTION  FACT(N) 

GIVES  FACTORIAL  OF  AN  INTEGER,  I.E.  N  •  (N-1)  «  (N-2)  • 
FACTORIAL  OF  ZERO  IS  1 

FACT= 1 

IF(N.LT.2)  RETURN 
Kzl 

DO  110  Is2,N 
K=K«I 
1 10  CONTINUE 
FACT=K 
RETURN 
END 
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FIG.  3  UNFILTERED  FLIGHT  DATA  -  (3)  LINEAR  VELOCITIES  AND  ACCELERATIONS, 

AND  MISCELLANEOUS  QUANTITIES 
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FIG.  4  NOISE  POWER  SPECTRA  -  (1)  CONTROL  INPUTS  AND 
DIRECTLY  LINKED  DISPLACEMENTS 
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FIG.  5  NOISE  POWER  SPECTRA  -  (2)  ANGULAR  DISPLACEMENTS  AND  VELOCITIES 
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FIG.  6  NOISE  POWER  SPECTRA  ~  (3)  LINEAR  VELOCITIES  AND  ACCELERATIONS, 

AND  MISCELLANEOUS  QUANTITIES 
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FIG.  7  EFFECT  OF  LOW  CUT-OFF  FREQUENCY  (=  1.12  Hz)  FILTER  ON 
UNFILTERED  FLIGHT  DATA 
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FIG.  8  EFFECT  JF  LOW  AND  HIGH  CUT-OFF  FREQUENCY  FILTERS  ON  VERTICAL 
ACCELERATION  SIGNAL  IN  REGION  OF  A  SUDDEN  CHANGE  ‘ 
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